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Abstract. We consider Schrodinger operators on metric cones whose cross section is a 
closed Riemannian manifold (Y, h) of dimension d — 1 > 2. Thus the metric on the cone 
M = (0, oo)r X y is dr^ + r^h. Let A be the Friedrichs Laplacian on M and Vb be a smooth 
function on y, such that Ay + Vo + (d - 2)^/4 is a strictly positive operator on L^iY), 
with lowest eigenvalue /iq and second lowest eigenvalue /^f, with /io,Mi > 0. The operator 
we consider is-ff = A + Vo/r^, a Schrodinger operator with inverse square potential on A/; 
^-H notice that H is homogeneous of degree —2. 

We study the Riesz transform T = V H~^^^ and determine the precise range of p for 
04 which T is bounded on L^[M). This is achieved by making a precise analysis of the operator 

(— I (H + 1)^^ and determining the complete asymptotics of its integral kernel. We prove that if 

|— { V IS, not identically zero, then the range of p for boundedness is 

d d \ 

^min(l + f + A*o, d) ' max(| - ^o, 0) j ' 
while if V is identically zero, then the range is 

57 fl, d 

<^ \^ max(f-^i,0)^ 

(-H The result in the case V identically zero was first obtained in a paper by H.-Q. Li [34j . 

1. Introduction 

^ The Riesz transform T on the Euchdean space M is defined by 

^ T = ^\,L 

fSj where Ajgd is the Laplacian operator. In this paper we study the Riesz transform T in a more 

\^ general setting of metric cones. A metric cone M is of the form M = Y x (0,oo), where 

(y, h) is a compact Riemannian manifold with dimension d—1. The cone M is equipped with 
Cs| the conic metric g = dr"^ + r^h. The Euclidean space M provides the simplest example of 

. . a metric cone, with cross section Y = S'^"^ with its standard metric. General metric cones 

enjoy a dilation symmetry analogous to that of Euclidean space, but no other symmetries in 
^ general. 

H The Laplacian on the cone expressed in polar coordinates is 

(1) ^ = -dl-—dr + \^Y, 



2000 Mathematics Subject Classification. 42B37, 58J05, 35J05. 

Key words and phrases, metric cone, Schrodinger operator, Riesz transform, inverse square potential, 
resolvent. 

The first author was supported by an ARC Future Fellowship FT0990895 and an ARC Discovery Grant 
DP1095448. 

1 



2 



ANDREW HASSELL AND PEIJIE LIN 



where Ay is the Laplacian on the compact Riemannian manifold Y. Then the Riesz transform 
T on the cone M is defined by T = VA^a^ where V is shorthand for (9r,r~^Vy), or in 
other words we measure the gradient on the cone using the metric g. The question of the 
boundedness of the Riesz transform on cones, i.e. for what p the operator T is bounded on 
LP{M), was answered by H.-Q. Li in |33j- The characterisation of the boundedness, stated in 
Theorem is in terms of the second smahest eigenvalue of an operator involving Ay . We 
provide a different proof to this result in Section [5] of this paper. 

Theorem 1.1. Let d > 3, and M be a metric cone with dimension d and cross section Y. 
The Riesz transform T = VA~2 is bounded on IJ'{M) if and only if p is in the interval 

max(^ -;Ui,0))' 

where > is the square root of the second smallest eigenvalue of the operator Ay + {^^)^ ■ 

More significantly, the methods used in this paper to prove Theorem |1.1| can be applied 
to study the boundedness properties of a more generalised class of operators, obtained by 
introducing an inverse square potential to the Laplacian. Let Vq :Y ^ Che a function on Y 
satisfying the condition 

(3) Ay + Vb(y) + (^)2>0 
and define 

(4) H = A + y^. 

Notice that H is homogeneous of degree —2, like the Laplacian. Condition ([s]) ensures that 
H is a strictly positive operator, so H~^/^ is well-defined. We can then define the Riesz 
transform T of the Schrodinger operator H by 



(5) 



Notice that ([3| allows our potential ^ = ^ to be "a bit negative" ; in particular, it allows Vq 
to be any constant greater than —(d — 2)^/4. 

The goal of this article is to find the exact interval for p on which the Riesz transform T 
with an inverse square potential V = Vo/r^ is bounded on IJ'{M), where M is a metric cone 
with dimension d> 2>. 



A necessary condition, stated in Theorem 1.2 for the boundedness was found in [28j by 



C. Guillarmou and the first author, in a slightly different setting — asymptotically conic 
manifolds. These are complete Riemannian manifolds {M°,g) such that M° is the interior of 
a compact manifold with boundary, M, which has a boundary defining function x for which 
the metric g has the form 

dx'^ h{x) 
x^ ' 

in a collar neighbourhood of dM, where h{x) is a family of metrics on dM. Here r = l/x 
behaves like the radial coordinate on the cone over dM: the metric in terms of r reads 
g = dr'^ + r^/i(l/r), so is asymptotic to the conic metric dr'^ + r^/i(0) as r — )• oo. In [28] , 
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potentials of the form V £ x'^C^{M) were considered; that is, the potentials decay as r ^ at 
infinity, and the limiting 'potential at infinity' Vq was defined by Vq := x~'^V\gj^. 

Theorem 1.2. i^ |28l Theorem 1-5]) Let d> 3, and {M°,g) be an asymptotically conic mani- 
fold with dimension d. Consider the operator P = + V with V G x^C°°(M) satisfying 

(6) AoM + Vo+ (^^) >0 on L'^iY), where Vq = ^ 

Let fiQ > be the square root of the lowest eigenvalue of the operator Suppose that P has 
no zero modes or zero resonance and that Vq ^ 0. Then VF~^^'^ is unbounded on LP{M) if 
p is outside the interval 

(7) ( ^ , ^ 

ymind + 1 + /Mo,d) max{^ - fj.o,0) 

The counter-example used in [28] to show the unboundedness of the Riesz transform can be 
easily adapted to the context of metric cones, so a similar result also holds for metric cones. 
Therefore the task now is to find a sufficient condition for boundedness. We will see that the 
sufficient condition involves the same interval ([T]) as in Theorem 1.2, so this interval gives us a 



complete characterisation of the boundedness of T with V ^ 0. Our main result is as follows. 



Theorem 1.3. Let d > 3, and M be a metric cone with dimension d and cross section Y . 
Let Vq be a function on Y that satisfies Ay + Vb(y) + (^^)^ > 0. The Riesz transform T 
with the inverse square potential V is bounded on L^{M) for p in the interval 

f d d \ 

(8) -j , , 

\^min(l + I + ;Uo,(i) max(| - /io, 0) y 

where /Uq > is the square root of the smallest eigenvalue of the operator Ay + Vo{y) + i^^)"^ . 

Moreover, for any V ^ 0, the interval ^ characterises the boundedness of T , ie T is 
bounded on LP{M) if and only if p is in the interval (^. 

Remark 1.4. If we specialize to positive potentials, i.e. V >0 and V ^0, then /xq > {d — 2)/2, 
and we see that the lower threshold for L^ boundedness is 1, and the upper threshold is always 
greater than d. On the other hand, for negative potentials V, i.e. y < and V ^ 0, the lower 
threshold for the boundedness is always greater than 1 and strictly less than 2, while the 
upper threshold is strictly less than d but strictly larger than 2. 

Remark 1.5. Part of these results are implied by a recent paper of Assaad [4j dealing with 
more general classes of potentials on M'^ or on complete Riemannian manifolds; see the end 
of Section 11.21 for further discussion. 



An immediate application of Theorem 1.3 is to show that the converse of the second part 



of \28\ Theorem 1.5], i.e. Theorem 1.2, is also true. As noted in |28| Remark 1.7], Theorem 



1.3 is exactly the missing ingredient. Therefore we have the following result. 



Theorem 1.6. Let {M°,g), P and /xq be as in Theorem | Then the Riesz transform 
is bounded on L'P{M) if and only if p is in the interval ([T]). 



A special case of Theorem |1.3| is the following result on the Riesz transforms with constant 
non-zero Vq, in which the boundedness interval is written in terms of the constant. 
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Corollary 1.7. Let M be a metric cone with dimension d > 3 and cross section Y. The 



Riesz transform T = V(A + 
and only if p is in the interval 



where c > 



'd-2\2 



(9) 



1.1. 



(10) 



2d 



and c 7^ 0, is hounded on L^{M) if 



2d 



ymin(d + 2 + yf{d - 2)2 + 4c, 2d) max(d - ^'{d - 2f + Ac, 0) j 

Strategy of the proof. Using functional calculus, we get the following expression, 

„ 2 



vr 



ViH + \^)-^d\. 



Because of the homogeneity of H, we obtain the resolvent kernel for {H+}?')~^ from {H+l)~^ 
by scaling the variables. So it suffices to analyze P := {H+l)~^ . We do this on a compactified 
and blown up space, which is designed so that the asymptotics of its kernel in different regimes 
can be understood. Let us use y as a local coordinate on the cross section Y . We particularly 
want to distinguish the diagonal behaviour of the kernel P~^{r,y,r' ,y'), from the behaviour 
as r or r' tend to zero or infinity. If we consider the kernel as living on [Y x [0,oo])^, as 
in Figure [T| then this has the defect that the diagonal meets the boundary hypersurfaces 
{r = 0}, {r' = 0}, {r = oo} and {r' = oo}, making the different asymptotic behaviours 
difficult to distinguish. To remedy this we perform blowups, as in [28j. As noted in that 
paper, the operator rPr is elliptic as a 6-differential operator near r = 0, that is, an elliptic 
combination of the '6-vector fields' rdr and dy.. On the other hand, as r — )• oo, P is an 
elliptic scattering differential operator, which is to say that it has an expression that looks 
like the Euclidean Laplacian in polar coordinates as r — )• oo, being an elliptic combination 
of dr and r~^cL^. Correspondingly we perform the 6-blowup (used to define the 6-calculus 
— see Section pi) for small r, that is, blow up the corner r = r' = 0, while for large r we 
perform two blowups (used to define the scattering calculus), namely we first blow up the 
corner r = r' = oo, followed by the boundary of the lifted diagonal at r = oo, obtaining the 
space illustrated in Figure [l} Now the diagonal is separated from the boundary hypersurfaces 
in Figure [l] and on this blown-up space, we can more easily construct the kernel of P~^ and 
describe the different types of asymptotics. 



r=r=M 




r=r'=0 




Figure 1. The space {Y x [0,oo])^, left, and the blown-up double space after 
three blowups, right. 
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Because the kernel behaves differently in different parts of the blown-up space, and es- 
pecially because we use different calculi near the two hypersurfaces zf and sf, we break the 
blown-up space into different regions, and construct the resolvent kernel in each region sepa- 
rately using different tools and techniques. In the end we patch up the constructions in these 
different regions to obtain the overall resolvent kernel. This construction of the resolvent 
kernel of H, ie the kernel of P~^, is done in Section [4| 

In Sectionjsj equipped with the knowledge on the behaviours of the kernel of at different 
parts of the blown-up space, we determine the boundedness properties of the Riesz transform 
T. Using a smooth partition of unity on the blown up space, we perform the integral ( |10[ ) 
and then break the kernel of T up into a near-diagonal part and an off-diagonal part. The 
near-diagonal part is a Calderon-Zygmund kernel and is bounded on for all p £ (1,00), 
while the off-diagonal part is bounded on typically a smaller range of p determined by the 
leading asymptotic behaviour at the boundary hypersurfaces marked 'Ibz' or 'rbz' in Figure[T] 



1.2. Relation to previous work. Cones have been studied since the 19th century, par- 
ticularly the problem of wave diffraction from a cone point which is important in applied 
mathematics, for example in [9j by A. Sommerfeld. Other notable early papers include [24| 
and by F. G. Friedlander and [10] by A. Blank and J. B. Keller. The Laplacians defined 
on cones were studied by J. Cheeger and M. Taylor in [16j and [17]. Many papers have been 
written about spaces with cone-like singularities. For example, the Laplacian and heat kernel 
on compact Riemannian manifolds with cone-like singularities has been studied in |18j by J. 
Cheeger and in [35] by E. Mooers, in [12J, J. Briining and R. Seeley studied the Laplcian on 
manifolds with an asymptotically conic singularity, and in [36j R. B. Melrose and J. Wunsch 
study the wave equation and diffraction on spaces with asymptotically conic singularities. 

The classical case of the Riesz transform on the Euclidean space M*^ goes back to the 1920s, 
and the case of one dimension (the Hilbert transform) was studied by M. Riesz in [40J. The 
paper [13] by R. S. Strichartz is the first paper that studies the Riesz transform on a complete 
Riemannian manifold. In T. Coulhon and X. T. Duong proved that the Riesz transform 
on a complete Riemannian manifold, satisfying the doubling condition and the diagonal bound 
on the heat kernel, is of weak type (1, 1), and hence is bounded on for 1 < p < 2. Since 
then, there have been many studies of the Riesz transform, of which we mention just a few: 
studies of the Riesz transform on complete Riemannian manifolds include [20], [7], [S], ^T\ : 
on Lie groups include [2], [3], [22], [Jl]; on second order elhptic operators [11], [23] . 

Many papers have been written on Schrodinger operators with an inverse square poten- 
tial. We only mention a few of the most relevant ones here. In f33], X. P. Wang studied 
the perturbations of such operators. In [14J, G. Carron studied Schrodinger operators with 
potentials that are homogeneous of degree —2 near infinity. In |13| by N. Burq, F. Planchon, 
J. G. Stalker and A. S. Tahvildar-Zadeh, the authors generalise the corresponding standard 
Strichartz estimates of the Schrodinger equation and the wave equation to the case in which 
an additional inverse square potential is present. In [31] the first author and A. Sikora in- 
vestigated one-dimensional Riesz transforms, including with inverse square potentials, with 
respect to measures of the form r'^~^dr, thus mimicking the measure on a d- dimensional cone. 

Now we turn to past results on the boundedness of the Riesz transform T with a potential 

V on metric cones. We have already mentioned the result (Theorem 1.1) of H.-Q. Li for 

V = 0, and the work [28] of C. Guillarmou and the first author on asymptotically conic 
manifolds. The method from [28] was based in part on the paper [15]. In [29] the two authors 
performed a similar analysis but allowed zero modes and zero resonances. In [6], P. Auscher 
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and B. Ben Ali obtained a result on M , stated in Theorem 1.8, which involves the reverse 
Holder condition. It is an improvement of the earlier results by Z.W. Shen in |42] . 

Theorem 1.8. Theorem 1.1]) Let 1 < q < oo. // V is in the reverse Holder class Bg 
then for some £ > depending only on V the Riesz transform with potential V is bounded on 
LP(M'^) forl<p<q + e. 

The reverse Holder condition V & Bq, implies that ^ > almost everywhere and V G 
^locO^'^)- ^ positive inverse square potential is in Bq if and only if q < d/2. So this theorem 



gives boundedness for 1 < p < d/2, which is smaller than the range obtained in Theorem 1.3 
for positive inverse square potentials (of course this is a very small subclass of i?q-potentials) . 

Very recently, Assaad and Assaad-Ouhabaz have proved results for Riesz transforms of 
Schrodinger operators which include some of our results. The following result is from [1]: 

Theorem 1.9. Let M be a non-compact complete Riemannian manifold with dimension 
d> 3. Suppose that the function y < satisfies A + (1 + £)V > 0, the Sobolev inequality 

(11) ll/ll^^(,,,^<l|V/||..(M), 

holds for all f € C^{M), and that M is of homogeneous type, ie for all x G M and r > 0, 

^,{B{x,2r)) <fi{B{x,r)), 

where fi is the measure on M. Then the Riesz transform T = V(A + V)^^ is bounded on 
L'P{M) for all p in the interval 



(12) 



2d 



d + 2 + {d- 



This result can be directly compared with ours in the case of Schrodinger operators of the 
form (A + where the constant c satisfies — (^^)^ < c < 0. In that case the lower threshold 
in ([9]) given by Corollary 1.7 is the same as the lower threshold in (12) given by J. Assaad's 



result. Also in it is shown that the Riesz transform for Schrodinger operators with po 
tentials in L'^/'^ on a d-dimensional Riemannian manifold obeying the Sobolev inequality (11) 



are bounded on (l,d) provided that this is true for the Riesz transform with zero potential. 
Note that this case just fails to cover inverse square potentials on cones, which are in 
Further results on signed potentials are proved by J. Assaad and E. Ouhabaz in 

2. Review of the 6-calculus and the scattering calculus 

In this section we briefly recall the key elements of the 6-calculus and the scattering calculus 
that we require in Section |4) For more details, see [37] or |27j for the 6-calculus, and |i38l [39] 
for the scattering calculus. 

2.1. 6-calculus. Let X be a manifold with boundary and with boundary defining function x 
(that is, dX = {x = 0} and 7^ at dX). The 6-calculus is a "microlocalization" of the set 
of 6-differential operators, namely those generated over C°^{X) by vector fields tangent to 
the boundary of X; near dX such vector fields are a linear combination of vector fields xdx 
and dy^, in terms of a local coordinate system (x, yi, . . . , Vd^i) with (yi, . . . , yd-i) restricting 
to a local coordinate system on dX. 

It is convenient to regard such operators as acting on 5-half densities, that is, multiples of a 
half-density taking the form \dx/xdyi . . .dyd-i\^^'^ near the boundary. Correspondingly, the 
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Schwartz kernels of such operators may be written as a distribution tensored with a 6-half 
density in each of the left and right variables. 

To define the 5-calculus, we first blow upj^X^ along (dX)'^ to obtain the blown-up manifold 
X'^ = [X'^; (dX)^], called the 6-double space. This produces a manifold with corners which 
has three boundary hypersurfaces: one defined by x/x' = (here and below we use the 
convention that unprimed variables on the double space are coordinates on the left copy, and 
primed variables are coordinates on the right copy), one defined by x' /x = and one defined 
by x + x' = 0. These are usually denoted lb, rb and ff respectively, but in accordance with our 
notational conventions for the space in Figure [l| we will call them Ibz, rbz and zf here (the 
'z' stands for 'zero' here and refers to the fact that the 6-blowup takes place at r = r' = 0). 

The small b-calculus ^'™(X), m G M, is defined as the set of 6-half-density-valued distribu- 
tions u on X| satisfying 

(i) u is conormal of order m with respect to diag;,, smoothly up to the hypersurface zf; 

(ii) u vanishes to infinite order at Ibz and rbz. 

Using the Schwartz kernel theorem, we interpret these as operators on (smooth functions) on 
X; the space "^^{X) extends to a bounded operator on L^. We also define 

m 

such operators are simply smooth 6-half-densities that vanish at Ibz and rbz. 

The 6-calculus is closed under composition; see [371 Prop 5.20] for the proof of the following 
result. 

Proposition 2.1. If X is a compact manifold with boundary then 

'i>]^iX) o ^f{X) C ^^+™'(X), 

where m, m! G M. 

Since our purpose is to invert elliptic 6-differential operators, it's important to know about 
parametrix constructions under the small 6-calculus. It is analogous to [33^ Theorem 18.1.24]. 

Proposition 2.2. If P is an elliptic partial differential operator of order k, then there exists 
an operator G in the small b-calculus of order —k such that 

Id-PG G ^j;°°(X), Id-GP G *^°°(X), 

and G with this property is unique up to an element of ■ 

For the proof, see \?>1\ Sec. 4.13]. This inversion property is not good enough for Fredholm 
theory, as the error terms Id —PG, Id —GP may not be compact. To investigate when an 
element in the small 5-calculus is compact, we introduce the indicial operator. 

Definition 2.3. Let A G ^'^(X) be a 6-pseudodifferential operator. The indicial operator 
hi^) is defined to be the restriction of the Schwartz kernel of A to zf. 



Here and below we use 'blow up' to mean real blow up; as a set, the manifold [X; S] obtained by blowing 
up X at the submanifold S is obtained by removing S and replacing it with its inward pointing spherical 
normal bundle. It is endowed with a difFerentiable structure that makes polar coordinates around 5* smooth 
functions on the blown up space. 
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The indicial operator /fc(^) can be interpreted as a translation-invariant operator on the 
cyhnder dX x M. As such it is an algebra homomorphism: 

h{PA) = h{P)Ib{A). 

The compactness of an operator is linked to its indicial operator. 

Proposition 2.4. Suppose that X is a manifold with corners, and A G with m < 0. 

Then A is compact on if and only if Ib{A) = 0. 

When inverting an elliptic partial differential operator in the small 6-calculus, the error 
term will usually have a non-zero indicial operator, and therefore will not be compact. In 
order to obtain an error term whose indicial operator vanishes, we have to expand the small 
6-calculus into a bigger calculus, called the full 6-calculus, in which the Schwartz kernels 
are permitted to have polyhomogeneous conormal expansions, i.e. expansions in powers and 
logarithms, at the boundary hypersurfaces Ibz, zf and rbz. 

To define polyhomogeneous cornormal functions, we need the notion of an index set. This 
is a discrete subset F C C x No such that every 'left segment' F ri{{z,p) : Re z < N}, N gR. 
is a finite set. Also, it is assumed that if {z,p) G F and p > q, q G N, we also have {z, q) E F. 

Given a boundary hypersurface and an index set, we can define polyhomogeneous conormal 
functions with respect to it. They are functions behaving like sums of products of powers and 
logarithms in one (and hence any) boundary defining function. 

Definition 2.5. Let X be a manifold with boundary and let H be its boundary. Given an 
index set F, a smooth function u defined on the interior X° of X is called polyhomogeneous 
conormal as it approaches the boundary H with respect to F if, on a tubular neighborhood 
[0, 1) X H of H, one has 

u{x,y) ^ az^p{y)x'' \og^ X 

{z,p)(iF 

as X — )• with a^^p smooth on H. Here, ~ means that the tail of the series, 

u = u- ^ az,p{y)x'' logP x, 

{z,p)eF,Rez<B 

is conormal and vanishes to order x^^"" for some e > 0, in the sense that |Vi . . . Viu'\ < Cx^^"^ 
for any finite number of vector fields Vi tangent to H applied to u' . 

Given a manifold with corners X, and an index family £ for it, i.e. an assignment of 
an index set for each boundary hypersurface, we define polyhomogeneous conormality of 
u G C°^{X°) by requiring that at each boundary hypersurface, u has an expansion with 
respect to the corresponding index set with coefficients that are polyhomogeneous conormal 
on the hypersurface; this sets up an inductive definition. See ^37^ Sec. 5.22] for details. 

Definition 2.6 (Full 6-calculus). The full b-calculus on X, where m is a real number and 
£ = (E'lbz) Erhz) is an index family for X^, is defined as follows. A distribution u on X^ is in 
^'™'^(X) if and only if u = ui + U2 + with 

(i) ui is in the small calculus 

(ii) U2 is polyhomogeneous conormal with respect to the index family (i?ibz, C°°, £^rbz), 
where C°° := {(n, 0) : n G No} is the C°° index set, and the index sets Sibz, C°° and 
E'rbz are assigned to the three boundary hypersurfaces Ibz, zf, rbz correspondingly; 
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(iii) U3 = /3*v, where /3 : X'^ — t- is the blow-down map and v is poly homogeneous 
conormal with respect to the index family £. 

Proposition 2.7 ([37, Prop. 5.46]). The full b-calculus on X is a two-sided module over the 
small b-calculus, i.e. 

and 

where m,m' £ M, and £ is an index family. 

The reason to introduce the full 6-calculus is that within it, we can construct parametrices 
of elliptic 6-differential operators with compact error term. For the proof of the following 
proposition, see [371 prop 5.59]. 

Proposition 2.8. Let P be an elliptic b- differential operator of order k whose indicial operator 
Ib{P) is invertible on L^{dX x M). Then there exists G in the full b-calculus of order —k such 
that the Schwartz kernels of the error terms E = Id —PG and E' = Id —GP are smooth across 
the diagonal, vanish at zf and are poly homogeneous conormal at Ibz, rbz with positive order 
of vanishing there. This implies that E,E' are compact on L'^{X). Necessarily (in view of 



Proposition 2.4), we have 

(13) h{G) = hiP)-'. 

2.2. Scattering calculus. Let X be a manifold with boundary dX and with local coordi- 
nates x,yi, ...,yd-i near dX, where x is a boundary defining function of dX. A smooth vector 
field V on X is a scattering vector field if it is x times a vector field on X, ie it has the form 

V = aox^dx + aixdy^ H h ad-ixdy^_^, 

with the coefficients oq, fld-i are smooth functions of x and y. Written in terms of r = x~^, 
these take the form 

V = -aodr + —dy, + ... —dy,_, . 

A scattering differential operator is one that is generated over C°°{X) by scattering vector 
fields. A key example is when X is the radial compactification of M°': then any constant 
coefficient vector field on M"^ is a scattering vector field viewed on X, and therefore any 
constant coefficient differential operator on M*^ is a scattering differential operator on X. The 
idea of the scattering calculus is to 'microlocalize' this set of differential operators. 

To define it we first need to blow up the product X"^ to produce the scattering double space. 
This is done in two stages: the first is to create the 6-double space X'^ = [X^; (dX)"^] as in the 
previous subsection. After this blowup, the diagonal lifts to be a product-type submanifold in 
X^, i.e. can be expressed as the vanishing of d coordinates in a coordinate system. The second 
step is to blow up the boundary of the lifted diagonal. The new boundary hypersurfaces so 
created are denoted bf and sf , respectively. 

Proposition 2.9. The interior of the scattering face sf in the scattering double space X|^ is 
a bundle over dX, and each fibre 0,y, y G dX , has a natural vector space structure. Moreover, 
any scattering vector field lifts from either the left or the right factor to be tangent to sf , and 
to be a constant coefficient vector field on each fibre. 
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It is convenient to regard elements of the scattering calculus (defined in the next para- 
graph) as acting on smooth scattering half-densities, i.e. taking the form at the boundary 
f\r'^~^drdy\^/^ , f G C°°{X). Thus the Schwartz kernels of such operators will be distribu- 
tions tensored with the half-density factor 

(14) r'^-^r"^'^drdydr'dy' . 

Definition 2.10 (Scattering calculus). The scattering calculus ^^''(X) of order (m, /) is de- 
fined as the set of distributions on X^^, times (14), satisfying 

(i) x~^v is conormal of order m with respect to the diagonal (more precisely the diagonal 
lifted to X"^^ uniformly up to sf, where x is a boundary defining function for sf; 

(ii) V vanishes to infinite order at the other boundary hypersurfaces. 
The order m is called the differential order of and / the boundary order. 

Remark 2.11. Using the Schwartz kernel theorem, elements of ^^''(X) may be interpreted 
as operators on half-densities on X. A scattering differential operator of order m acting on 
half-densities is in ^^c^{X). 

The scattering calculus is closed under composition. 

Proposition 2.12. |39l Eqn. 6.12] Let X he a manifold with boundary, and m,l,m',l' £ M, 
then 



Like Proposition 2.2 on the parametrix constructions under the small 6-calculus, under the 
scattering calculus we also have a result analogous to |33| Theorem 18.1.24]. 

Proposition 2.13. Suppose that P G '^sc{X) is elliptic. Then there exists G G ^^c'^{X) 
such that 

PG-Id, GP -lde^-^^^{X). 



Similarly to the case of the indicial operators in Section |2.1[ the normal operator of A G 
\E'"^'*'(X), denoted Nsc{A), is defined to be the restriction of the Schwartz kernel of A to 
the scattering face sf. This restriction can be interpreted (in a canonical way) as a smooth 
function on dX valued in densities on each fibre. These densities can be interpreted as 
convolution operators on functions (or half-densities) on each fibre. Under this interpretation, 
normal operators can be composed, and the action of taking normal operators is an algebra 
homomorphism: 

Proposition 2.14. [Ml Eqn. 5.14] Let A and B be elements o/^'*'°(X). Then 

Nsc{AB) = NsciA)Nsc{B). 

As with the indicial operator, vanishing of the normal operator is related to compactness: 

Proposition 2.15. Let A G ^''"'^(X) with m < 0. Then A is compact if and only if Nsc{A) 
vanishes identically. 

Remark 2.16. Alternatively, we may describe the boundary behaviour in the scattering calcu- 
lus by taking the fibrewise Fourier transform of each convolution operator, obtaining a family 
of multipliers; this is known as the normal or boundary symbol. Composition in terms of the 
boundary symbol is simply pointwise product. 



RIESZ TRANSFORM ON METRIC CONES 



11 



Proposition 2.17. If A G 'i'"^'^{X) is elliptic with invertible normal operator, then there 
exists B G such that E = AB - Id is in i.e. its Schwartz kernel is 

smooth across the diagonal and rapidly vanishing at the boundary of X^^. In particular, E is 
compact and hence A is Fredholm, with parametrix B. Necessarily, we have 



As discussed in the Introduction, we will construct the resolvent kernel = (i? + l)~^ on 
a compactified and blown up version of its natural domain M^, using both h- and scattering 
blowups. We start by compactifying in each factor separately, i.e. we pass to the compact 
space [0, oo]r x y x [0,00]^' x Y, where [0,oo] indicates the compactification of [0,oo) by a 
point at r = 00, such that 1/r is a boundary defining function at r = 00. As noted in the 
Introduction, rPr is an elliptic 6-differential operator down to r = 0, while P itself is an 
elliptic scattering differential operator up to r = 00. Therefore we perform the 6- blowup at 
r = r' = and the scattering blowups at r = r' = 00. This means that we blow up the 
corner r = r' = 0, the corner r = r' = co and finally the boundary of the lifted diagonal 
{r = r' ,y = y'} at r = r' = 00. 

We label the boundary hypersurfaces of [0, oo]r x y x [0, 00],.' x y by Ibz, Ibi, rbz and rbi 
according as they arise from {r = 0}, {r = 00}, {r' = 0}, or {r' = 00}, respectively. The 
new boundary hypersurfaces created by blowup are labelled zf, bf and sf, according as they 
arise from the blowup of r = r' = 0, r = r' = c« or the boundary of the lifted diagonal at 
r = r' = 00, respectively. The resulting space after the blow-ups at r = r' = and r = r' = 00 
is called the blown-up space. See Figure [l| 

We next discuss local coordinates near the various blown up faces. Near zf, local coor- 
dinates are {r/r' ,r' ,y,y') when r/r' < C (that is, away from rbz) and {r,r' /r,y,y') when 
r'/r < C (that is, away from Ibz). Near bf and away from sf the situation is similar: coordi- 
nates are (r'/r, r'~ ,y,y') for r'/r < C and {r/r',r~^,y,y') for r/r' < C. Near the interior of 
sf, coordinates are (r — r',r{y — y'),y,r^^). In the case that M is Euclidean space M*^, with 
Euclidean coordinate z, then z — z' is a. linear coordinate on each fibre of sf (cf. Proposi- 
tion [2^. In particular, the diagonal is defined hy r /r' = l,y = y' for small r (that is, away 
from sf ) and r — r' = 0, r{y — y') = or r — r' = 0, r'(y — y') = for large r (that is, away 
from zf). The following result about the diagonal will be useful later. 

Proposition 3.1. Let if : [0, 00) — )• [0, 1] be an increasing smooth function such that f{x) = x 
for X G [0, 2] o-nd (p{x) = 1 for x G [1, 00). Then the function 



where z = (r, y) and z' = (r', y'), is a quadratic defining function for the diagonal in the blown- 
up space,' that is, Odiag 

> 0, the diagonal lifted to the blown up space is given by {odiag = 0}, 
and the Hessian of Odiag directions normal to the diagonal is positive definite. 

Proof. The formula for the distance on a metric cone is given by 



NscB = {NscA) 



1 



Proof. See [38l Section 6]. 



□ 



3. The blown-up double space 




(15) 
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(The second line is because when dY{y,y') > tt the fastest way to get from (r, y) to {r',y') is 
to go straight to the cone point and back out again.) So near the diagonal we have 

d{z,z'f = (r - r'f + 2rr' (l - cos ((iy(y,y'))) 
(16) / ^ 

= (r - r'f + rr' Uy (y, y'f + 0(dy (y, y'f)\. 

Near the sf-face, we have 

adiag(^, z') = d{z, z'f = {r- r'f + rr' (^dy (y, y'f + 0(dy (y, y'f)^ , 

which is good for a quadratic defining function for the diagonal. To see that we recall from 
the discussion before this proposition that near sf the diagonal is defined by r' — r = 0, and 
r{y — y') = or r'{y — y') = 0, and we also recall the standard fact that dY{y,y'f is a 
quadratic defining function for the diagonal of Y'^ for any closed Riemannian manifold Y. 
Near the zf-face, we have 

adiag(., z') = = - 1)' + ^ {My, y'f + 0{dy{y, y'f)^ , 

which is again good for a quadratic defining function for the diagonal, as here the diagonal is 
instead defined hy ^ = I and y = y' . □ 



3.1. Densities on the blown- up space. By a smooth 6-half-density on the blown- up space 
we mean a half-density of the form 



u{r, r',y,y'] 



dr dr' . 

i-dydy 

r r 



where u is smooth. (This is perhaps misleading since it is only a 6-half density in the usual 
sense away from sf . However, we shall only use this when either r or r' is small, in which case 
it certainly is a 6- half density.) Let x = ^,x' = ^. Then by a smooth scattering- half-density 
we mean a density of the form. 



v{x,x',y,y'] 



dxdx'dydy' 



where v is smooth. In terms r and r' it becomes, 



.1. 1 



v{r, r', y, y')|r'^+V"^+^d(-)d(^)dydy' 



, 1 

/I 2 



v{r, r', y, y')\r'^-^r"^-'^drdr' dydy' 



, 1 

/I 2 



The scattering half-density \r r' drdr'dydy'\^ is a bounded nonzero multiple of the 
Riemannian half-density. We will usually consider the resolvent as acting on Riemannian 
half-densities, in which case the kernel of itself is a Riemannian (distributional) half- 
density on the blown-up space. However, when wc study the properties of a kernel near the 
zf-face, we write it as a 6-half-density; this is more natural in view of the fact that we use the 
6-calculus near zf . 
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4. Resolvent Construction 

4.1. The operator H. Let M be the metric cone over (Y, h). The Laplacian on the cone M 
expressed in polar coordinates is 

where Ay is the Laplacian on Y. This operator is positive and symmetric on the domain 
C^{Y X (0, oo)), i.e. smooth functions supported away from the cone tip. The operator A is 
defined to be the Friedrichs extension of this symmetric operator. 
For any function Vq : Y ^ C, we define the operator 

This is a natural class of operators: as both A and are homogeneous of degree —2, the 
operator Hvq has the same homogeneity. For simplicity of notation, we write Hvq simply as 
H. The following proposition tells us for which Vq is the operator H positive. 

Proposition 4.1. Suppose that Ay + Vo{y) + (^^)^ is a positive operator on L'^(Y). Then 

the operator H is also positive. 

Proof. Wc work in polar coordinates, and consider the isometry U : L^{M;r'^~^drdy) 
Lp'{M;r^^drdy) defined by 

(17) Uf = r-2f. 

Now for / G L'^{M;r~^drdy), we compute 

UHU-^f = ri(-dl- '-^dr + lAy + y^yif. 



fci)^vi(,))J,/ + l*/-a,V + ^AK/. 



Therefore we have, 



(^j(-(r5r)^ + Ay+(^)Vyo(y))j)/ 
= - dr(rdr{\f)) + ^Ay/ + [[^)" + V,{y)^ 1/ 

=dr{\f) - dlf + lAy/ + ^C^y + Voiy)) 1/ 

= - ^/ + l^rf - dlf + lAy/ + ((^)' + Fo(y)) ^/ 

'^'-'^^o(y)U/ + ^a./-a.V + ^Ay/ 



=UHU-^f^ 
We have established 



(18) 



= 1 ( - {rdrf + Ay + Vo{y) + (^^^)') l 
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Make a substitution s = Inr, then the space L'^{M;r ^drdy) becomes L'^{M] dsdy), and we 
have 

?7ifC/-i = e-^ f - + Ay + Vb(y) + " 



2 

From here we can clearly see that the operator H is positive if Ay + Vo(y) + (^^)^ > 0. This 
completes the proof. □ 

Remark 4.2. Note that as we have d > 3, the condition Ay + Vo{y) + (^^)^ > means that 
the potential y = ^ is allowed to be "a bit negative" . 

4.2. The Riesz transform T. Suppose we have a function Vq on y which satisfies Ay + 
^oiy) + (^^)^ > 0. Then the Riesz transform T with the inverse square potential V = ^ is 
defined to be 

where the size of the derivatives are measured using the cone metric g, i.e. derivatives of 
bounded length are given by {dr,r~^dy-). 

Our aim is to find out the precise range of p for which the Riesz transform T is bounded 
on U'{M). Following [15] and |28j. we do this using a 'resolvent approach' as opposed to 
the more common 'heat kernel approach'. Using functional calculus, we have the following 
expression, 

9 /-oo 

T=- V{H + \^)-^d\. 
71" Jo 

We see from this equation that in order to understand T, we need to know the properties of 
{H + A^)^^. Because H is homogeneous of degree —2, we only need to compute {H + 1)~^, 
then use scaling. Let P = H + 1; we will proceed to study P~^ . 

4.3. A formula for the resolvent. We now proceed to find an explicit formula for P~^. 
However as we will discuss later, the formula has good convergence properties in only certain 
regions of the blown-up space. From Equation ( 18 ) we have 



P = H + l = - {rdrf + Ay + Vb(y) + + f V^" 



Let P' denote the differential operator consisting of the terms in the middle. That is, 

2 . . . 9 /d- 2\2 



(19) P' = -{rdrY + Ay + Vo{y) + r'+ ^ ^ 

We take P' to act on half-densities, using the flat connection that annihilates the Riemannian 
half-density \r'^~^drdh\^^'^ on M. Now let P be the differential operator given by the same 
expression (19), but endowed with the flat connection on half-densities annihilating the 6-half 
density \dr/rdh\^^'^. Since U maps this 6-half density to the Riemannian half-density, these 
two differential operators are related by 

(20) P = U'^P'U. 
Therefore, 

(21) p = r~^Pr~^. 



Since P is self-adjoint. Equation ( 21 ) shows that P is also self-adjoint. (Note that for operators 
on half-densities there is an invariant notion of self-adjointness, since the inner product on 
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half-densities is invariantly defined.) Denote G = P^^, G = P^^; the Schwartz kernels of G 
and G are related by 

(22) G = rr'G. 

(Again, we emphasize that this is an identity involving half-densities: if we write G = 
K\{rr'Y-^drdr'dhdh'\^l'^ and G = k\{rr')'^drdr' dhdh'\^l^ then we have 

(23) K = {rr'f-'^'^K.) 



So we just need to determine G, then Equation (22) gives us G. 



We now proceed to work out an expression for G. Let {fj.'jjUj) be the eigenvalues and 
corresponding L^-normalized eigenfunctions of the positive operator Ay -|- Vo{y) + (^^)^- 
We also let Ilj denote the projection onto the u^-eigenspace. Then we have 

(24) p = ^n,r„ 

j 

and 



r 
J 

where 

(25) fj = -{rdrf + + ^| = -r'^d^ - rdr + 

As in |28j , the kernel of the inverse of Tj is written in terms of modified Bessel functions 
Ifij{r) and Kfj_.{r) (see [H Sec. 9.6]) in the form 

' I^^{r)K^^ir')\fd^\K r<r', 
^M,(r)/^, 



fr\r,r') 
We know that 



K,Ar)I.Ar')\^^\', r>r', 



j 

hence in terms of the kernels, we have 



(26) G(r r' y y') = [ iy)u, jy')!,, {r)K,^ (rQ | f '^dhdh' | ^ , r < r' , 

\E,MyMy')K,Ar)lHi^TT'^dhdh'\K r>r', 

where dh denotes the Riemannian density with respect to the metric on Y. While this is an 
exact expression for G, it is not a very useful expression near the diagonal, as it has poor 
convergence properties. Therefore we shall glue it together with a pseudodifferential-type 
parametrix in order to determine its properties close to the diagonal. However, sufficiently 
far from the diagonal, the series has very good convergence. We proceed to show this. 



4.4. Convergence of the formula. By the symmetry of (26), it suffices to consider the 
region {r < r'}; here we work with the sum, 

(27) j;u,(y)^/^^(r)i^^^,(r'). 
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Prom [U Sec. 9.6], we have representations 



7r2 9^^r''^ /"°° , 1 



and 



r(/^ + 5) A 

We now estimate each of these integrals in a way that is uniform as /i — t- oo. When r < 1, 
the integral in the expression for is uniformly bounded in /x > 0, and hence we see that 

(28) \W\ < + 1/2) "^^^"^ - ^' 

where C is independent of r and /i. On the other hand, for r > 1, we estimate e~''r^/^/^(r): 

TT2T{fl+^) J-1 

- r(/x + i) 

< c V^-^ / e - 

- r(/x + i) io r 

< c 

" r(/x + i) 

with C independent of /u, for /i > 1/2. This gives rise to an estimate of the form 

(29) \W\ < ^f(^^7^Y/^ when r > 1. 
We next estimate -fC^ in a similar way. For r < 1, we estimate 

nh-^r^ /-oo _^ ^ 



r(/i+ 2) Ji 



(30) 



r(/^ + ^)yo 



< C ^ \ / e-H^^-^r-^^' dt 



^2-^r-^r(2/u) 

r(// + ^) 
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On the other hand, for r > 1, we compute 



r(/i+2) 



1 

„ln-n^H+l/2 roc 

' e-\t{2r + t)Y-2r-^^^dt 



where we made a substitution t — )• r[t — 1) in the integral. We now estimate 
which gives rise to an estimate 



(31) 



K,{r)\<C^p^m^^[{2rY-^'^T{^, + l/2),T{2^,)) 



Now to absorb the factor in the first argument of the maximum function, we sacrifice 

half of our exponential decay: we estimate e~''/^r^~^/^ by bounding it by the value where it 
achieves its maximum in r, which is when r = 2/i — 1: 

g-r/2^M-l/2 < e-(2A'-l)/2(2^ - 1)^-1/2 < ^6-^2 V^^^-^/^) < C2^'T{^l). 



Then we can use this in (|31|) to estimate 
(32) 



< C 



-r/2 -^t 



r(/^ + 5) 

Finally using the identity (see [1, 6.1.18]) 
(33) 

we obtain 



max 



(22A'r(/x)r(^ + i/2),r(2/i) 



22/i-l 1 

r(2/x) = ^r(^)r(/. + J 



(34) 



< Ce-'^''^r-^'2^^'T{^). 



Hence, when r' > Ar, I^{r)K^{r') is bounded above by 
(35) 

We emphasize that the constant C is independent of /x > 1/2, r and r' here. Noting that the 
combination of F factors is uniformly bounded in each case (using (33) again), we find that 
for r' > 4r, I^{r)K^{r') is bounded above by 



|-^(.). 2-^^r(2,) o<r<r'<l 
(r(M+i/2)) 



(36) 



'C{^)\ 0<r<r'<l 
C7(^)^e--72, r<\<r' 

c(^Y^-'''\ 1 < r < r'. 
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By Hormander's L°°-estimate, see [32], we know that ||uj||oo ^ ^/^j^ • Therefore each 

term in the series is bounded above by C ^'^^^ (jr)^^ ^ . To continue the discussion on 
convergence, we need the following lemma. 

Lemma 4.3. Suppose that are the eigenvalues of Ay + Vo(y) + (^^)^, then for any 
< /3 < 1, and any M, N > 0, the sum 

^J,j>M 

converges for all < a < P, and it is bounded uniformly in a. 
Proof. Note that for any iij > 2M, we have 

/.,-M = M + (^,-2M)>M+(^^l^) = |. 

Therefore, 

tMj>2M Hj>2M fij>2M 

There is an integer A^i(/3, N) > 2M such that for all j > Ni{/3, N), < /3~4. It follows that 

E s E "f + E r 

tMj>2M 2M<^ij<Ni(l3,N) H>^iW'^) 

(37) <\{j:fi,<N,{P,N)}\Ni{P,Nf+ ^ 

tij>Ni(/S,N) 

fij>Ni{l3,N) 

where the constant C > comes from the Weyl's estimate, which states that for any fj, > 1, 
we have 

(38) |{i:/ii<M| <CA'- 

We continue to estimate the part of summation greater than Ni(f3,N). An implication of 
(38) is, for any j £ N, we have 



Therefore, 



1 

J \ d-l 
C 
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There is N2{f3,C) G N such that for aU j > iV2(/3,C), we have H^)'^-'^ > log^j, where 
7 = > 1. Then 



4VC'' < 



(39) 



0<i<Af2(/3,C) i>Af2(/3,C) 

<iV2(/3,C)+ 

i>7V2(/3,C) 

<iV2(/3,C) + ^. 
The remaining part of the summation is from M to 2M, 

Y ii^jOi^''^'^ < \{j ■■ /Wj < 2Af}|(2M)^ < Ci2M)'^-'^{2M f = C{2Mf+^-\ 

M</ij<2A/ 

Bringing all the parts together, we have 

2 

(40) Yl /^f < C(2M)'^+^-i + CNi{{3, iV)'^+^-i + N2{f3, C) + < oo. 

Note the finite constant depends on M, N, C, f3 but not a, therefore we have uniform bound- 
edness in q. □ 




Figure 2. Support of (5/ 



Proposition 4.4. The expansion (27) is polyhomogeneous conormal at Ibz. 



Proof. Since the functions Ifj,{r) and K^{r) have expansions in powers at r = (including 
logarithms in the case of when ^ is an integer), the individual terms in the series are 
polyhomogeneous conormal. So consider the tail of the series. Lemma |4.3| implies that the 
sum of the tail of the series, that is over //j > M is bounded by Cr*^e~^^ for small r. We 
can apply the same argument to derivatives of the series. In fact, the derivatives of 1^ and 
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can be treated as above, showing that for r < 1, fj, ^{rdr)^I^ and /i ^{rdr)^K^, and for 
r > I, fi~''d^Ifj_ and fi'^'d^K^ satisfy the same estimates as and K^. Moreover, we have a 
Hormander estimate ||V*-'^'''Uj ||oo < CkfJ^i^ i)/2+fc derivatives of Uk- Thus derivatives only 



give us extra powers of fi, which are harmless as Lemma 4.3 applies with arbitrary powers of 



II. □ 



Proposition 4.4 implies, in particular, that G decays exponentially, with all its derivatives, 
as r' — >■ DO, ie when approaching the boundary rbi. Similarly in the region p- > 4, as r — ?■ oo, 
ie when approaching Ibi, the kernel is also exponentially decreasing. Therefore we cut off G 
to restrict it away from the r = r' to obtain a well defined operator Gj with the kernel 

(41) Gf{r,r',y,y') = G{r,r' ,y,y')(^x[^) ^'^(v) 

Here x is a smooth cutoff function x '■ [0, oo) — )■ [0, 1] such that x([0, 1/2]) = 1 and x([l) c«)) = 
0. Thus the support of G/ is contained in {r/r' < 1/4} U {r/r' > 4}, as illustrated in Figure 
[2| (The subscript 'f stands for 'far from the diagonal'.) 
At last, similar to (22), we define 

(42) Gf = rr'Gf. 

4.5. Near Diagonal. The formula obtained in the previous section doesn't converge near 
the diagonal, so in this section we construct an operator Gnd, which is good near the diagonal. 
The subscript nd means "near diagonal" . 

Near the zf-face we consider the 6-elliptic operator P. In order to keep it away from the 
sf-face, we multiply it with a cutoff function, so we consider Px{r), where x '■ [0, oo) — )• [0, 1] is 
a sm ooth cutoff function as above. By the ellipticity of P near the zf-face, and by Proposition 



2.8 



there is G^^^ in the full 6-calculus such that 



(43) PG'Jax{r) = xir) + E,f, 

where E^f is smooth across the diagonal and vanishes to first order at zf (as a 6- half density). 
Let 

(44) G^f, = rr'&J,, 
then we have 

PGl[x{r) = x{r) + E,f 

where E^f = {r' /r)Ezf is smooth across the diagonal and vanishes to first order at zf as a 
6-half density. 

Near the sf-face the operator P is elliptic in the scattering calculus. We multiply it with 
1 — x{f) to keep it away from the zf-face, ie we consider the operator P(l — xii"))- Since 



^'(l - X{r)) 
Proposition 



is elliptic near the sf-face, and its normal operator A^n + 1 is invertible, by 
.17, there is G^'^ in the scattering calculus such that 

PG'i^,{l-x{r)) = l-x{r) + Esf, 

where the error term Egf is smooth across the diagonal and vanishes to infinite order at the 
boundary of the blown up space. 
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Now we define Gnd by 

,4r^ Ar' , 
Then we have 

PGnd = Id +End, 

where the error term End is smooth across the diagonal, and vanishes to first order at zf (as 
a 6-half density) and to infinite order at all other boundary hypersurfaces. We may assume 
that Gnd is supported close to the union of the diagonal, zf , and sf . 
We now define our global parametrix to be 

(45) Ga = Gf + Gnd. 

4.6. The indicial operator at zf. In this subsection we show that the leading behaviour 
of Ga at zf agrees with that of Gnd- To do this, it suffices to show that the indicial operator 
of G agrees (at least for r/r' < 1/4 and r/r ' > 4, where we have shown convergence of the 



2.8 



series) with that of Gnd- Proposition 
Let us now determine this indicial operator 
The indicial operator of P is 



the indicial operator of G^^ is equal to Ib{P) ^- 



h{Px{r)) = -{rdrf + Ay + V^iy) + 



Let fj,j,Uj,Ilj be the same as defined in Section 4.3 Here, instead of (24) and (25) we have 



where 

Similar to Section 



4.3 



Si 



then the kernel S, is 



h{Px{r))=Y,U,S,, 

-{rdr 



Sj ^(r,r') 



2^ij V r' ) 
2iJ,j V r ) 



r r 

dr dr' 
r r' 



r < r', 
r > r' . 



Hence 
(46) 



[h{P))-\s,y,y') 



dhdh 



1 

'I 2 



s > 1, 
s < 1, 



The convergence of this sum can be analyzed using Lemma 4.3 the sum converges smoothly 
for s < 1 and for s > 1. 

Now we determine the leading behaviour of G at zf . We only consider the case ^ <1 /4, 
as the case > 4 is completely parallel. Recall from expression (26), for ^ < 1/4 we have, 

G{r,r',y,y') = ^Uj{y)uj{y')I^^{r)K^^{r') 



dr dr' 



dhdh' 



We use the limiting forms for small arguments from |lj Sec. 9.6], that is when r, r' — )• 0, 
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and 
(48) 

Therefore, 

(49) Go{r,r',y,y') 



-1 



- + U[r I logr I 



r \H 
.1 



dr dr' 



r r 



dhdh' 



+ 0(r'^ I log/ 1 



Since 



dr dr 



ds dr 



(46) of {Ib{P)) ^{s,y,y') = hiG^,) when restricted to the z/-face. 



expression (49) of G{r,r' ,y,y') at zf is indeed consistent with expression 



Finahy, since the cutoff function used to define Gj is x(4r/r') +x(4r'/r), and that used to 
define Gnd is 1 — x(4?'/r') — x(4r'/?^), we see that Ga has the same leading asymptotic at zf 
as G^-^, namely h{Py^ ■ 



Gf + G 



nd: 



4.7. Construction of P ^. We have constructed an approximate inverse Ga 
let E be the corresponding error term: 

PGa = ld+E. 

We next try to solve away E to obtain our final G = P^^. We begin by summarising the 
properties of Ga and E. 

Proposition 4.5. As a multiple of the Riemannian half-density \r'^~^r'^^^drdr'dhdh'\2 on 
the hlown-up space, the kernel Ga is the sum of two terms. One is Gnd, supported where 
1/8 < r/r' < 8, and is such that p'^^'^Gnd is conormal of order —2 with respect to the 
diagonal uniformly up to both zf and sf , where p^f is any boundary defining function for zf , 
and is rapidly decreasing at bf. The other term G f = Ga — Gnd satisfies: 

(i) it is smooth at the diagonal, and polyhomogeneous conormal at all boundary hyper- 
surfaces; 

(ii) it vanishes to infinite order at Ibi, rbi and bf ; 



(iii) it vanishes to order 1 

(iv) it vanishes to order 2 



2 + Po at Ibz and rbz; 
d at zf . 



Proof. The properties of Gnd follow from properties of the full 6-calculus and of the scattering 
calculus recalled in Section [21 

The diagonal part of property (i) of Gf is clear: in fact it is supported away from the 
diagonal. Polyhomogeneity of Gf at Ibz and rbz follows from Proposition 4.4 and the sym- 
metry of Gf, while polyhomogeneity (in a trivial sense, with an empty index set) at Ibi, rbi, 
and bf follows from the exponential decrease of Gf as r or r' tend to infinity, as shown by 
Lemma 14.31 



We obtain the vanishing order at Ibz from equations (41 ) and (42 ). Since r is the boundary 



defining function for Ibz, we need to work out its power. Clearly one power of r comes from 
(42), while Ifigir) in (41) gives us the power r^o. Then the difference between the 6-half 
density and the Riemannian half-density gives us a power of r~'^/^ (as in (23)). Combining 



these we conclude that the vanishing order at Ibz is 1 — ^ + /^o- The vanishing order at rbz 
is similar. 

Last, we show (iv). Since both r and r' vanish at zf, to obtain the vanishing order of Ga at 



zf, as a scattering-half-density, we combine the powers of r and r in (44) and (42) with the 
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factor (rr') '^1'^ involved in the change from a 6-half density to the Riemannian half-density. 

d 
2 



So the order of vanishing isl — 5 + I — 1 = 2 — d. □ 



Proposition 4.6. The error term E has the following properties on the blown-up space: 

(i) it is smooth in the interior; 

(ii) it vanishes to the first order (as a h -half- density, or to order 1 — d as a Riemannian 
half- density) at the zf-face; 

(iii) it vanishes to infinite order at Ibz, rbz, Ibi, rbi, sf and bf; 

(iv) it is compact on L?'{M); in fact its Schwartz kernel is Hilbert- Schmidt. 

Moreover, the k-fold composition satisfies similar conditions, with (ii) strengthened to 
vanishing to order k at zf as a h -half- density. 

Proof. Property (i) follows from the choice of Gnd- Property (ii) follows from the fact that 
the indicial operator of {rr')~^Ga is equal to Ib{P)~^, as shown in the previous subsection. 
Property (iii) follows from the fact that elements of the scattering calculus vanish to infinite 
order at bf, the fact that Gj is equal to the exact inverse of P outside the region {1/8 < 
r/r' < 8} (so in fact E is supported in this region, hence vanishes in a neighbourhood of Ibz, 



rbz, Ibi and rbi), and the exponential vanishing of Gf as either r — )• cx) or r' — )• 00 — see (36). 
Properties (i), (ii) and (iii) show that E has an kernel, proving Property (iv). 

To show the last remark, we use a smooth cutoff function to divide E into two parts, 
E = Eb + Esc, where Eb is an order —00 operator in the 6-calculus, vanishing to first order at 
zf , and Esc is an order (—00, 00) operator in the scattering calculus. Then E^ = {Eh + Esc)''- 
Any mixed terms will vanish to infinite order at each boundary hypersurface. Of the remaining 
terms, using the composition properties of the b- and scattering calculus recalled in Section [2j 
E^ is order —00 in the 6-calculus and vanishes to order k at zf , while Eg^ is order (—00, 00) in 
the scattering calculus. Moreover, E^ is supported where {8^^ < r/r' < 8'^}, hence vanishes 
in a neighbourhood of Ibz, rbz, Ibi and rbi. □ 

We proceed to solve away E. To achieve that, we would like to invert Id+S. But it might 
not be invertible: if not, we perturb Ga so that ld-\-E becomes invertible. 



Since E is compact on L'^{M), according to Proposition 4.6, ld-\-E is Fredholm of index 
0, and its null space and the complement of its range both have the same finite dimension, 
say A^. Removing the null space gives us an invertible operator, and to achieve that we add 
a rank operator to Ga- To construct the rank A^ operator we need the following lemma. 

Lemma 4.7. There exist smooth functions ipi, ...,ij)N, <pN on M such that 

(i) ipi, ...jipN span the null space of Id -\-E , and P(j)i, P(f)]\f span a space supplementary 
to the range of Id -\-E; 

(ii) they are 0{r°°) as r — t- and 0{r^°°) as r ^ 00. 

Proof. We choose the tpi to be any basis of the null space of Id -\-E. To obtain property (ii) 
for the ipi, we note that = hence iterating, we have t/'j = E'^^ipi for each A^ > 1. 

Now we consider mappi ng p roperties of the operator E^ . First, writing E = E^ -\- Egc as in 



the proof of Proposition 4.6, it is easy to see that Esc and VEsc map L'^{M) to {r)~^ Lp'{M) 
for arbitrary K. (Here V is shorthand for the vector of derivatives {dr,r~^dy.).) As for E^, 
since it has negative order in the 6-calculus and vanishes to first order at zf, we see that Eb 
maps L?'{M ) to rL?'{M). Since the kernel (r/r' YE has the same properties as E listed in 



Proposition 4.6, it follows that Ei, maps r"'L?'{M) to r'^~''^L^(M) for any a. Also, applying a 
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derivative V = {dr,r~^dy.) to E^, it is still of negative order in the 6-calculus, though no longer 
vanishing at zf, so we see that \7Eb maps r°'L'^{M) to r°'L'^{M) for any a. Summarizing, we 
have 

E boundedly maps r^L^iM) r"+^(r)"-^L^(M), 
(50) „ J. ^ 

VE boundedly maps r^L^iM) r^ir)-^ L^{M). 

Applying these properties of E iteratively, we see that E'^^ maps L?'{M) to {r)~'^^ (M) 
for any A^. Hence, using Sobolev embeddings, tp is smooth and has rapid decay both as r — )• 
and r — )• oo. 

As for the (pi, to show that we can choose functions (pi, . . . , (pN as above, it is sufficient to 
show that the range of P on the subspace S of smooth half-densities satisfying (ii) is dense 
on L?'{M). If this were not true, then there would be a nonzero half-density / G L?'{M) 
orthogonal to the range of P on such half-densities: that is, we would have 

{Pu, f) = 0, for ah n G S. 

Since S is a dense subspace, this implies that Pf = distributionally. By elliptic regularity 
this means that / is smooth and Pf = strongly, but since P is invertible on this implies 
/ = 0, a contradiction. Therefore we can choose the (pi £ § as desired. □ 

Let Q be the rank A^ operator 

N 

Q = ^(Pi{'^i^ •)> 

i=l 

where {ipi, ■) means the inner product with ipi. The functions ipi, (pi, ■■■4'N are chosen 

as in Lemma Then we have 



P{Ga + Q) = ld+E + PQ, 
which is invertible. From here we obtain 

P-^ = {Ga + Q)(Id +E + PQ)-\ 



Using property (ii) of Lemma 4.7, we see that Ga + Q has the 'same' properties as Ga, ie 



it has those properties listed in Proposition 4.5, and E' := E + PQ has properties (i) - (iv) 



listed in Proposition 4.6 Define operator S by 

S= {Id+E'Y^ -U. 

Then we can write 

P-^ = {Ga + Q){ld+S). 
We need to know the properties of S. 



Lemma 4.8. The operator S has properties (i) - (iv) listed in Proposition 4-6 
Remark 4.9. A similar analysis was done in |30| Sec. 5.4]. 
Proof. Using the identities (Id +S')(Id +^') = (ld+E'){ld+S) = Id, we obtain 
(51) S = -E' + E'^ + E'SE'. 
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For any positive integer N, we substitute the expression (51) into itself 2N — 1 times, and we 
get 

AN 

(52) S = ^{-lyE'^ + E'^^SE'^^. 

Using the last part of Proposition 4.6, we see that the term X]j=i(~l)"'-^'"' prop- 
erties listed in the Lemma, so we focus on the term Sn := E'"^^ SE''^^ . Using (50), we see 



that r ^(r)^^vi^'*i?'^^ and r' ^ {r')'^^\/[^^ E'"^^ are bounded operators on L^. Since S 

is Hilbert-Schmidt, it follows that {rr')"^ {r)^^ {r')^^vi^^v[^^ Sn has an kernel. Using 
Sobolev embeddings, this gives regularity and vanishing (at the boundary) of Sn of some 
finite order N + 0(1), and hence the same finite order regularity and vanishing of S. Since 
this argument can be made for any A^, this proves that S has the properties (i) - (iv) listed 



in Proposition 4.6 □ 



To summarise, we have 

G = P-^ = {Ga + Q){ld+S), 
where Ga + Q has those properties listed in Proposition |4.5[ Id +5 is a compact operator. 



and 5 has those properties listed in Lemma 4.8 Our final step is to analyze the composition 



{Ga + Q)(id+S) and show that G itself satisfies all the conditions listed in Proposition 4.5 
We summarise key information about G = obtained through our construction in the 



following theorem. To state it, define u = 1 — ;\;(4r/r') — x(4r'/r) where x is as in (41 ); thus, 
a; is a smooth function on the blown up space supported away from Ibz, Ibi, rbz and rbi, and 
equal to 1 on a neighbourhood of the diagonal. Also let pz{ be a boundary defining function 
for zf. 



Theorem 4.10. Let Gc = ujG and = (1 — oj)G. Then, as a multiple of the Riemann- 
ian half- density, ie the scattering-half- density \r'^^^r"^^^drdr'dhdh'\2 , on the blown-up space, 
p'^Y'^Gc is conormal of order —2 with respect to the diagonal uniformly up to both zf and si. 



while Gs satisfies properties (i)-(iv) listed in Proposition 4-5 



Remark 4.11. The subscripts c and s are chosen to indicate that Gc is the part of G which is 
conormal at the diagonal, while Gs is the part of G which is smooth at the diagonal. 



Proof. We have already proved these properties for Ga, in Proposition 4.5, so we need to 
check them for the terms Q + QS + GaS = G — Ga- Since Q and QS both are smooth and 
vanish to infinite order at the boundaries, these terms trivially satisfy all the conditions. So 
it remains to check that GaS has the same properties as Ga- 

We write GaS as a sum of two parts. Let r] : [0, oo) — )• [0, 1] be a smooth cutoff function 
such that 7/([0, 1]) = and rj{[2,oo)) = 1. The first part r]{r)Ga'f]{r') is in the scattering 
calculus. Note that r]{2r')S is also in the scattering calculus, and that 

{7]{r)GaV{r')){v{2r')S) = {r,{r)GaV{r')) S- 



Therefore by Proposition |2. 12] this term is in the scattering calculus. The second part Ga — 
r]{r)Ga'f]{r') is in the full 6-calculus. (Although the support of this term meets the boundary 
hypersurfaces Ibi and rbi, its Schwartz kernel is rapidly vanishing there, enabling us to regard 
it as living in the 6-calculus.) In a similar sense, S is in the small 6-calculus (it vanishes rapidly 



at every boundary hypersurface except zf). Therefore by Proposition 2.7, {Ga—r]{r)Gar]{r')) S 
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is in the full 6-calculus, with the same index sets at Ibz and rbz as Ga- Therefore the required 
properties for Gc and property (i) for Gs follow. Also, since S vanishes to first order at zf, 
the same is true for the composition (Ga — r]{r)Gar]{^'))S. So GaS has the same vanishing 
orders (or better) at the boundary hypersurfaces as Ga- □ 



Remark 4.12. In the case of the potential y = 0, we have /Uq 



2 — 1. So the vanishing 
order in item (iii) of Theorem 4.10 becomes 0. It's consistent with the case when the cone is 
R"^ and the potential V = 0, when the cone tip can be chosen arbitrarily, and G is smooth 
everywhere. 



The vanishing orders of G 
in Figure [3j 



P ^ at various boundaries of the blown-up space are shown 




Figure 3. The vanishing properties at various boundaries 



Remark 4.13. The construction of G = in this section is sketched in the paper [28] by 
C. Guillarmou and the first author, but details are lacking. It's not fully justified in [28j that 
the kernel is in the scattering calculus near sf, and in the 6-calculus near zf. For this reason 
we have given complete details in this section. 

5. The Boundedness of the Riesz transform 

5.1. Estimate on the kernel. Recall that the Riesz transform T with the inverse square 
potential V = defined in Section |4| can be expressed as 

T = - V{H + x^y^dx, 

^ Jo 

where H is given by Q, and recall that H is homogenous of degree —2. Our analysis of the 
Riesz transform will be based on the following estimate on the kernel T{z, z'). 

Proposition 5.1. We have the following estimate on the kernel ofT, 

POO 

\T{z,z')\< / A'^-2|v(G(Az,Az'))|^^A, 
Jo 

where G = P^^ = {H + 1)^"*^, with properties listed in Theorem 4-10. 
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Proof. This comes from the relationship between (H + A^)~^ and G = {H + 1)^^, which is 

{H + A2)-1(z, z') = X^-^{H + l)-^{\z, Xz'). 

The power —2 of A appears because H is homogenous of degree —2. Remember these kernels 
are Riemannian half-densities, and this accounts for the power d of A: 

\{XrY-^{Xr'f-^d[Xr)d{Xr')dhdh'\^ = X\'^~\"^-^drdr'dhdh'\L 

□ 

5.2. Boundedness on L'^{M). 

Proposition 5.2. The Riesz transform T with the inverse square potential V = ^ is bounded 
onL^{M). 

Proof Our assumption is Ay + Vb(y) + (^)^ > 0, ie A + ^Vb(y) > 0. Hence there is e > 
such that A + (Yz^Vb(y) > 0. That is, A + ^Vb(y) > eA. From here, 

(T/,r/) = (A(A + ^Vo{y)yh\ (A + V) 

< (£"'(A + ^Vo{y)) (A + ^Voiy))-h, (A + ^^0(2/))"^/) 

= ^"'((^ + ^My))h, (A + ^Voiy))-h) 

= e~HfJ). 

Therefore T is bounded on L2(M). □ 

5.3. The diagonal region. To understand the Riesz transform on L^, we break up G as in 



Theorem 4.10 Here we will write Gi for Gc = ojG (recall a; = 1 — x(4r/r') — x(4r'/r)), and 
we further decompose Gg = G2 + G3, where G2 = Gx(4r/r') and G3 = Gx(4r'/r). Notice 
that G2 and G3 are supported away from the diagonal, in particular where the infinite series 



(26) has good convergence properties as shown in Proposition 4.4 We correspondingly break 



up the Riesz transform into three pieces. Thus we have 



(53) 



Ti{z,z') = - / X''-^V,{Gi{Xz,Xz')) dX. 
Jo 



We now show that Ti is of weak type (1,1). For that we first need to estimate the derivatives 
of Gi. 

Lemma 5.3. Let d{z, z') denote the distance between z and z' on M . On the support of lo, 
we have pbf ^ d{z., z')~'^ , where /3bf is a boundary defining function for hi. 



Proof. Let z = {r,y) and z' = {r',y'). Observe from (15) that d{z,z') is bounded above by 



r + r' . Therefore in the region {1/8 < r/r' < 8} we have 

d{z, z')-^ > (r + r')-^ = r'-\l + > 
As r'~^ is a boundary defining function for bf on the support of lo, the result follows. □ 
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Lemma 5.4. The kernel Gi satisfies the estimate that for any integer j > 0, we have 

' d{z,z'f-^-i, d{z,z')<l, 
diz,z')-^, d{z,z')>l, 



\K^'G,{z,z')\< 
for any N > 0. 



Proof Note that Gi is supported in the region Ri = {1/8 < r/r' < 8}. Since p'^^'^Gi is 
conormal of order —2 with respect to the diagonal, by Proposition |3.1[ near the diagonal we 
have 

\pti'G,{z,z')\<a^^. 
Near zf, adiag = so it follows that near zf, we have 



p: 



zf 



\Gi{z,z')\<diz,z') 



l\2-d 



Away from zf, adiag = d{z,z')'^, therefore 



lag 

2-d 
2 

diag 



Gi{z,z')\<a^?=d{z,z'f-'. 



Now let's consider the behaviour of Gi near bf. By Theorem 4.10[ we know that it vanishes 
to infinite order at bf, while by Lemma 5.3, we know that /3bf ^ d{z,z')~'^. Therefore near 
the bf-face we know that 

\G^{z,z')\<d{z,zy, 

for any > 0. The rest of i?i is easy because after we take away the neighbourhoods near 
zf, bf and the diagonal, we are left with a compact set, on which both Gi and are 
continuous with d{z,z')~^ being non-zero. Therefore we can conclude that 



\Gi{z,z')\< 



d{z,z'f-^, d{z,z')<l, 
d{z,z')-^, d{z,z')>l, 



for any > 0. Then using the conormality of G at the diagonal and polyhomogeneous 
conormality of G at the boundary hypersurfaces, we obtain the estimates on V''^ ^,Gi{z, z'). 
□ 



Proposition 5.5. The operator Ti maps L^{M) into L}'^'^°-^{M) . 

Proof. We just apply the Calderon-Zygmund theory, see [261 Section 8.1.1]. It is sufficient to 
verify the following conditions: 

(i) Ti is bounded on L?{M); 

(ii) \T^{z,z')\<-^^,- 

(iii) \V,Ti{z,z')\ < ^-j+T and \V ,>Ti{z, z')\ < ^ 



{d{z,z')) {d{z,z')) 

for some constant C > 0. 



d+l 1 



We already know from Proposition 5.2 that T is bounded on L (M). So to verify condition 



(i), we just need to show T — Ti is bounded on L'^{M), which is covered by Proposition 
in Section [531 



5.14 
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Now we show conditions (ii) and (iii). By Lemma 5.4 we know the kernel Gi satisfies, with 
any A > 0, 



and 



\V,{Gi{Xz,Xz'))\ < X\{V,Gi){Xz,Xz')\ 



\vl{Gi{Xz,Xz')) \ < X^\{ViGi){Xz,Xz') 



< 



< 



' X'^-'^d{z,z'y~'^, Xd{z,z')<l, 

X-^+^d{z,z')-^, Xd{z,z')>l, 

\ X'^-'^d{z,z')-'^, Xd{z,z')<l, 

\X-^+^diz,z')-^, Xd{z,z')>l, 



for any > 0. We use this to estimate Ti{z, z') 



T,{z,z')\< / X''-^\v4Gi{Xz,Xz'))\dX 



< 



d(z,z') 



d{z,z'Y~'^dX+ I X''-''-'d{z,z')-''dX 



d-N-lj(^ J\-N , 



d{z,z')-'^ + d{z,z') 



d(z,z') 
oo 



x'^-^'-'dx 



d(z,z') 



d{z,z')-'^ + d{z,z') 



/\-d 



(Choose N = d+1. 



d{z,zY' 

Now estimate the derivative with respect to z. The z' case is similar. 

2 



|v,ri(z,z') 



TT 



< 



< 



X'^'^V ,V ,{Gi{Xz,Xz'))dX 
X'^^'^\vl{Gi{Xz,Xz'))\dX 

'^'^ ' d{z, z'^dX + / A'^-^dlz, zT^'dX 



d{z,z')-'^-^ + d{z,z') 



/\-N 



d{z,z') 
OO 



d(z,z') 



d{z,z')-'^-^ +d{z,z')~^-^ 
2 



(Choose N = d + 2.) 



d{z, z' 



This completes the proof. 



□ 



By interpolation, we obtain the following proposition. 

Proposition 5.6. The operator Ti is bounded on L'p(M) for any p > 1. 

Proof. By Marcinkiewicz Interpolation Theorem, we know that Ti is bounded on L^{M) for 
all 1 < p < 2. The same holds for the adjoint of Ti. Using duality, we get boundedness for 
1 < p < oo. □ 
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5.4. OfF-diagonal region. To study the boundedness of the two off-diagonal operators T2 
and T3, the following lemmas will be useful. They are similar to |311 Lemma 5.4] but not 
covered by it. 



Lemma 5.7. Consider the kernel K{r, r') defined by 

K{r,r') -- 

Ifa + /3 = d, (3>0, and p satisfies 



0, r > r' . 



(54) p < ^— rr, 

max(a, Uj 

then K is hounded as an operator on L^{'R+;r'^~^dr). 

Proof. The proof is essentially taken from [31j. To find out for what p the operator with 
kernel K{r, r') is bounded on L^'(R_|_, r'^~^dr), we consider the isometry M : LP{M.^,r'^~^dr) — >■ 
LP{R+,r-'^dr) defined by 

{Mf){r)=rif{r). 

Then the kernel of the operator K = MKM'^ : LP{R+, r^'^dr) LP{R+,r-'^dr) is 

K{r,r)=rpr pK{r,r) = { — ) ^X{-r<r'}- 
Perform a substitution s = Inr, s' = Inr', then K{s, s') is an operator on LP(R, ds), and 

J\[s, s ) — e p X{s-s'<o}- 

This is a convolution operator, so it is bounded provided the kernel is an L^-function with 
variable s — s'. Since s — s' < 0, we want —a + ^ > 0. That is, 

d 

p < 



max(a, 0) 

Since we want p > 1, we require a < d, ie /3 > 0. □ 



Corollary 5.8. Let K{r, r', y, y') he a kernel on the cone M satisfying 

\K{r,r\y,y')\ < 
If a + (3 = d, /3 > 0, and p satisfies 



r "'r' ^ , r < r', 
0, r > r'. 



(55) p< ^ 



max(a, 0) ' 

then K is hounded as an operator on U'{M]r'^^^drdh). 



Proof. It follows from Lemma 5.7 and the fact that the cross section Y has finite volume. □ 
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K{r, r') 



Lemma 5.9. Consider the kernel K{r, r') defined by 

Jo, r < r', 

If a + (3 = d, a > 0, and p satisfies 

(56) p> . f 

then K is hounded as an operator on LP(K^;r'^^^dr). 



Proof. By duality and Lemma 5.7 □ 



As before, Lemma |5 .91 has a corollary about the boundedness of operators on the cone M. 
Corollary 5.10. Let K{r,r' ,y,y') be a kernel on the cone M satisfying 

{n r < r' 

-a -/3 : ; 

If a + 13 = d, a > 0, and p satisfies 
(57) p> 



min(a, d) 

then K is hounded as an operator on U'{M]r'^^^drdh). 
Proposition 5.11. The operator T2 is hounded on LP{M) for 

P < 7^ ^' 

max(2 — A^O) Uj 

where fiQ > is the square root of the smallest eigenvalue of the operator Ay + VQ{y) + {^^)'^ . 



Proof We will use the following three boundary defining functions 

Pzf = r', pihz = Prbi = {r')~'^. 
r 

By Theorem 4.10[ we have 

(59) |G2(r,r',y,y')l < pIi" p\J^'' pfu = ,i-|+Mo/i-i-.o^/)-oo^ 

where ^0 > is the square root of the smallest eigenvalue of the operator Ay + Vo (?/) + ( ^^)^. 
From here we know that 

J';y2-c(j,l-f+/iO^/l-f-Ato A<- 
I. ' — r' ' 

for all > 0. That means, by polyhomogeneous conormality of G2, if /xq 7^ | — 1, we have 



(60) \V,{G2{Xr,Xr',y,y'))\ 
for all iV > 0. 



< 



2 
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Using (53) and the fact that G2 is supported where r < r', we estimate 



\T2{r,r\y,y')\< / A^-^j V, (G2(Ar, Ar', y, y')) |dA 



< 



Jo / 1 



1 



,N-di-N 



< ^-f +M0/-f-W for JV > ^0 + 



r \ Mo- 2 



By Corollary 5.8, we conclude that T2 is bounded on LP{M) provided that 

d 

p < 



□ 



Remark 5.12. When V = 0, then /xq 



1, and its first eigenfunction uq is a constant 



function. In Section 5.5 we will improve estimate (60) to obtain a bigger range for p for this 
special case. 

Proposition 5.13. The operator is bounded on LP{M) for 

d 

p > 



min(l + I +/xo,d) 

where fiQ > is the square root of the smallest eigenvalue of the operator Ay + VQ{y) + {^^)'^ . 
Proof We will use the following three boundary defining functions 

r' 



Pzf = r, Prhz 



-, Plhi = {r) ^ 



By Theorem 4.10 we have 



(61) 



\G^{r,r',y,y')\ < p.V'/'iV^'VS = r'-l'>^^V'-l+^^ {r) 



where ^0 > is the square root of the smallest eigenvalue of the operator Ay + Vo(y) + (^^)^. 



It follows that, as in the proof of Proposition 5.11, by the poly homogeneous conormality of 
G3, 



\V,{Gs{Xr,Xr',y,y'))\ < , 
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for all > 0. Then using ([53]) and the fact that G3 is supported where r' < r, we have 

/>oo 

mry,y,y')\< / A^-2|v4G3(Ar,Ar',y,y'))|dA 

JO 

1 1 



J. 2 t^'jf'^ 2 



-l-f-Mo/l-f+W) 



dA 



+ 



Applying Corollary 5.10 we conclude that T3 is bounded on U'{M) provided that 

d 

p > 



min(l + I + fj,o,d) 



□ 



Proposition 5.14. The operator Ti is bounded on L'^{M). 



Proo f. Since 2 satisfies the boundedness criteria in both Proposition 5.11 and Proposition 
5.13 the operator T2 + ?3 = T — Ti is bounded on Lp'{M). The operator T is bounded on 
L'^(M) by Proposition 5.2, and from here the boundedness of Ti on L'^{M) follows. □ 



Remark 5.15. Proposition 5.14| completes the missing part in the proof of Proposition 5.5 
5.5. Proofs of main results. 



Proof of Theorem \1.3[ Since T = Ti + T2 + T3, we just combine Proposition 5.6, Proposition 
5.11| and Proposition |5.13 to prove the first part of this theorem. 

For the second part, with V ^ 0, for p outside the interval ([s]), the counterexamples from 
|28^ Section 5.2] serve to show the lack of boundedness of T on LP{M). (For purposes of 
comparison, note that the variables x and x' in [28] correspond to ^ and p- in this paper.) □ 

Proof of Theorem \l.l\ Suppose that the potential V is identically zero; we proceed to show 
that the upper threshold for boundedness is p = d{d/2 — fii)^^- Notice that Ti and T3 
are automatically bounded on this extra range, so we only have to consider T2, which has an 
expression of the form 



(62) 



T, 



2 r°° 

{z,z') = -x(4r/r') / A^'^V, (G(Az, Az'))^^A. 



We recall that G = rr'G and substitute the infinite series (26) for G here, and consider the 



first term in this sum separately from the rest. Since no = d/2 — 1 when Vq = 0, the first term 
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here is (as a multiple of the Riemannian half-density — recall this gives us an extra factor of 
^rr')-'^/^, as in ^) 

(rr')i-'^/2no(y)uo(y')4/2-i(r)i^d/2-i(r-') 
When Vb = 0, the eigenfunction uo{y) is constant. Also, Id/2-iii^) = cr'^/'^^^ + 0{r'^^'^~^^) and 
is conormal at r = 0, implying that Vr-(r^~'^/^/^/2-i(^)) = 0{r). Hence 

V.(ri-'^/\o(y)/d/2-i(r)) = 0(r); 

that is, in this special case, applying the derivative Vz makes the kernel vanish to an additional 
order, instead of one order less as is usually the case. Therefore, after taking the gradient in 
the left variables, this term is bounded by 

r' > 1 

for any integer A^. Now we put this in (62 ) and find that the contribution to T2 of the /XQ-term 
is bounded by 




l/r' 



< Crr' 



l/r' 

-l-d 



X^-^lx^-'^rr'- 



f-N\ 
) 



dX 



5.8 



that this 



Remembering that this term is supported in {r < r'}, we see from Corollary 
term is bounded on LP for all p G (1, 00). 

So consider the remainder of the series. The argument in the previous subsection applies, 
except that the series now begins with the /ii term rather than the term, so we have 
boundedness in the range (58) with /xi replacing completing the proof. □ 
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